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Abstract 

The paper is devoted to the study of stabihty of equihbria of a 
delay differential equation that models leukemia. The equation was 
previously studied in [5] and [6], where the emphasis is put on the 
numerical study of periodic solutions. Some stability results for the 
equilibria are also presented in these works, but they are incomplete 
and contain some errors. Our work aims to complete and to bring 
corrections to those results. Both Lyapunov first order approximation 
method and second Lyapunov method are used. 
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1 Introduction 

The study of the mathematical model of periodic chronic myelogenous 
leukemia considered in [5], [6] may be reduced to that of the equation 



i(t) 



1 + x{t) 



where /3q, n, 6, k, r are positive parameters and k = 2e~"'^ , with 7 > 0. 
We do not insist here on the significance of the function x{.) or in that of 
the parameters, since these are extensively presented in [5], [6]. We only 
remind that the unknown function, x(-), should be nonnegative, being a 
non-dimensional density of cells. 

The two studies cited above are mainly devoted to the numerical inves- 
tigation of the delay differential equation above. The stability of equilibria 
study in [5], [6], that is reduced to a few lines, is, unhappily, incomplete and 
contains some errors. 
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Our work aims to correct the errors in the stabihty conditions presented 
in [5], [6] and to present some aspects concerning the dynamics generated 
by equation ([1]), aspects that were not pointed out there. 

In Subsection 1.1 we prove that the Cauchy problem associated to our 
equation has an unique defined on [— r, oo) bounded solution. In Subsection 
1.2, following [5], [6], the two equilibrium solutions, as well as the linearized 
equation and the characteristic equation for each of these, are presented. 
Section 2 deals with the stability study of the equilibrium points. We use 
the results of [5] and, in a first stage, we obtain results valid for both equi- 
libria. In the subsequent subsections we analyze the stability of the two 
points individually. It is important to perform this separate study since 
the conclusions are very specific to each equilibrium point. In Section 3 we 
comment the stability results in [5], [6], pointing out the errors therein. 

1.1 Existence and uniqueness of solution 

We make the notation B = C([— r, 0],M) (the space of continuous, real 
valued functions defined on [— r, 0], with the supremum norm, denoted by 
|x|o). Given a function x : [— r, T) i->- R, T > and a < t < T, we define 
the function xt € S by Xj(s) = x{t + s). 

Equation ([TJ may be written as 

x = h{xt), (2) 
where h : B M, and we impose to this equation the initial condition 

xo = (l)eB. (3) 

Remark that if the initial condition is a positive function, then x{t) 
can not become strictly negative. Indeed, let ti be the first moment when 
x{ti) = 0, (that is x{t) > for t < h). Then x{ti) = k^^0^^^ > 0, hence 
x{t) > for t > ti in a neighborhood of ti. 

We study the existence, uniqueness and domain of existence of solutions. 

The function h is globally Lipschitz. Indeed, by denoting 
f3{x) = /3o/(l + x"), we have 

!^(/3(x)x)| </3o(n + l), 
and thus, for any ipi, ip2 & B 

\hi^i)-h{ip2)\ < {(3oin+l)+6)\MO)-MO)\+k/3oin+l)\ipii-r)-M-r)\ < 

< [{k + l)Po{n + I) + 6]\ipi - ip2\o. 
From here the continuity of h follows also. 
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Theorem 2.3 of [3] implies that problem ([T]), ^ has an unique solution 
defined on an interval [0, T). 

If we take (^2 = 0, by using the fact that h{0) = 0, we obtain that 

\h{^i)\ < [(A: + l)/3o(n + l)+5]|v9iIo 

thus the function h is also completely continuous. 

We prove that the solution is bounded. For this, we multiply equation 
([1]) by x{t) and we get 

• / N / X ^ c- 9 / N , ^ x(t — r)x(t) 
x(t)x(t) = -/3o- - (t) + kj3o--^ / ^ \ < 

e- 2/ ^ ^ ^^(*) 1 , n x^(t-r) 



and 



2 l + x"(t-r) 2e l + x'^(t-r) 



dt e 
We chose an e > such that rj := 25 — ek/3Q > and we obtain by integration 

hence the solution of problem ([I]), (l3|) is bounded. Theorem 3.2 of [3] implies 
that the solution is defined on the whole positive real time semiaxis. 

Hence, for any cj) £ B, problem ([1]), ([3]) has an unique defined on M+ 
bounded solution. We can thus associate to this problem the semigroup of 
operators on B, {T{t)}t>o, T{t){(j)) = xt{(p), where x{t,(j)) is the solution of 
eq. ([1]) with initial condition xq = (j). 



1.2 Equilibrium solutions 

In this subsection we, inevitably, follow [5]. 
The equilibrium points of the problem are 

XI = 0, x2 = (^(fc-l)-l)^/^ 

The second one is acceptable from the biological point of view if and only if 
it is strictly positive that is, if and only if 

^(fc-l)-l>0. (4) 
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In terms of r, by using k = 2e "''^ , the above inequality may be written 

as 

r < Tmax ■■= In - ( 1 + — I , (5) 

and since the delay r is positive, the condition (5//3o < 1 follows. 

The biological interpretation of function /3 [5] shows that the condition 
/3(a^2) = S/{k — 1) > should be fulfilled. This is equivalent to A; > 1. 

The linearized equation around one of the equilibrium points is 

z{t) = -[B + 6]z{t) + kBz{t-r), (6) 

with B = (3'{x*)x* + x* = xi or x* = X2. 

The eigenvalues of the infinitesimal generator of the semigroup of op- 
erators generated by equation Q are the solutions of the characteristic 
equation 

X + S + B = kBe-^\ (7) 



2 Stability of equilibrium points 

In order to investigate the stability of the equilibrium solutions, we first 
try to establish the conditions in which all the eigenvalues have strictly neg- 
ative real part, in order to use the linear approximation Lyapunov method. 

We rely on the work that exhaustively solves the problem of finding 
necessary and sufficient conditions on the parameters such that the equation 
A = ai -|- 026"^ has only solutions with strictly negative real part. 

We denote 5 + B = kB = hence ([7]) becomes 

\+p = qe-^'-. (8) 

By taking A = ^ + iuj, and by equating the real, resp. the imaginary parts 
in our equation, we obtain 

fi + p = qe''^^ cos{ujr), (9) 

ui = —qe~^^ sin(a;r). 
It is useful to consider the case ^ = in the above equations, 

p = q cos{ujr), (10) 

ui = —qsm{ujr). 
The results in [2] imply the following 

Proposition All solutions A of eq. ^ satisfy ReX < 0, if and only if 
a) q < 0, < —pr < 1 and -p < —q < (wq 
or 
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b) q < 0, p > and -q < {ujq 
or 

c) q > 0, p > 0, and q < p, 

where uq is the solution in {0,-K/r) of the equation 

Ldcot{ujr) = —p. (11) 

Remark. If we divide the first equality in (|10|) to the second one, we 
obtain ([TT]) . Hence relations (fTO|) are equivalent to the set of relations (fTT]l 
and + = . 

In order to express cjq in a more direct form, we consider the function 
T : [0, vr) i-^ (-co, 1], given by 

ny) = { ll^o'""^' (12) 

The function is a bijection and we can equivalently define ojq, the solution 
of (mD, as 

coo = -T-\-pr). (13) 
r 

We express the conditions in Proposition in terms of r. We first remark 
that 

2^ 2 2 

cot (ti;orj cos'^ [ujQr) 

The two inequalities in of the above Proposition may be written as 

\P\ 



\P\ < \Q\ < 



cos(a;o?^)| 



Since p < 0, the solution ujq of equation ([TT]) is such that uor G (0, it/2). 
Hence the above inequality is equivalent to 

< - < 1, cos(a;o?^) < 
Q Q 

and the second one is equivalent to arccos(|) < ujQr < tt/2. To conclude, 
case a) is described by the inequalities 

arccos(^) i 

q<p<0, —<r<—. (14) 

^0 \p\ 



In case b) q < 0, p > 0, and we must have 

P 



-q < 



cos(a;or)| 
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In this case, a;ocota;or = —p < 0, and since ujQr £ (0, vr), we must have 
ujQr G (7r/2,7r). The above inequahty is equivalent to 

I cos{ujQr)\ < 

1^1 

and this one is satisfied if 

P 

p/\q\ > 1 or {p/lgl < 1 and (— cos(a;or) < — )}. 

-q 

The last condition is equivalent to 

p p TT p 

— 1 < - < Oand cos(a;or) > - ^ — < ujQr < arccos(-) 44> 
q g 2 q 

TT arccos(2) 

^ — < r < ^. 

2^0 ^0 

Remark. The case q < 0, p = corresponds to ujQr = 7r/2, and the 
eigenvalues lie to the left of the vertical axis if and only if —qr < tt/2. 

We can now translate the discussion above to our concrete problem. 
I. If B < 0, then two situations may occur. 

A. 6 + B < 0. In this situation, ReX < for all eigenvalues A if and only 
if \5 + B\ < \kB\ and 

arccos ( (6 + B)/kB) 1 , , 

where 

oJo = -T-H-iS + B)r), 
r 

with T given by ()12p . 

If the studied equilibrium point is X2, the condition r < rmax must be 
also fulfilled. 

B. 5 + i? > 0. In this situation, ReX < for all eigenvalues A if and only 

if 

6 + B>\kB\ori5 + B< \kB\ and r < =((i±^)«| (ig) 



where, again 



ujo = -T-H-{S + B)r), 
r 



with T given by 

II. If S > 0, then we can only have 6 + B > 0, and in this situation 
ReX < for all eigenvalues A if and only if 

kB <6 + B. 

Even if the above discussion seems comprehensive, it is still useful to 
consider the two equilibrium points separately and to discuss their stability. 
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2.1 Stability properties of xi 

In this case, B = (3q > 0, hence the necessary and sufficient condition 
for the negativity of the real part of ah eigenvalues is 



Since the condition ^ for the existence of the second equilibrium point, X2, 
is the reverse of the above inequality, it follows that xi is stable as long as it 
is the single equilibrium point. When the second equilibrium point occurs, 
xi becomes unstable. 

We inspect the eigenvalues at ^(/c — 1) = 1. Equation ([7]) in this case is 

X + 6 + l3o = kpoe-^\ 

and, since kf^Q = 5 + Pq, admits the solution A = 0. Hence the change of 
stability occurs by traversing the eigenvalue A = 0. 

2.1.1 Stability of xi when ^{k - 1) = 1 

In this case, the "first order approximation" theorem is of no use, since 
is the eigenvalue with greatest real part. We use a Lyapunov function in 
order to prove stability of the zero solution. 

However, since for our problem x{t) > 0, the concept of stability should 
be interpreted in the following way: 

for every e > there is a 6 > such that if </>(s) > 0, s G [— r, 0], and 
\(j)\Q < 5, then < x(t,(j)) < e for any t > 0, where, as above, x(t,(j)) is the 
solution of (CP with condition 

If y : ;B i-> M is continuous, the derivative along the solution x {■,(/)) of 
the Cauchy problem ([2]), ([3]) is defined as [3] 

y ((/>) = hm sup - y ((/.)]. 

Definition V is a Lyapunov function on G C B if V is continuous on 
G andV <0 on G. 

Theorem [1]. If V : B M is a Lyapunov function and there exist a 
continuous increasing function a : [0, oo) i— )■ [0, oo), with a(0) = and 

a{\m\)<V{cl>), 

then the solution x = of equation (0j is stable and every solution is 
bounded. 

We construct below a Lyapunov function for our problem, for the con- 
sidered parameter values. 
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Let us consider the function G{u) = f^2s/{l + s^)ds. We define 



rds. 



We have 



. (l + (/)"(s))2 

62(0) (/.2(-r) 



and by using the equality 

x(o, </.) = -mm) + mo) + mmrm-r), 

we obtain 



(l + (/."(0))2 1 + 0"(O) ' "■^'^"(l + (/."(0))(l + 0"(-r)) 



+ 



(l + (/)"(0))2 (l + 0«(-r))2 
from where, with the inequahty 



2<P{0)H-r) 



< 



'(0) 



+ 



</,2(-r) 



(1 + </'"(0))(l + </>"(-r)) - (1 + 0"(O))2 (1 + ./.'^(-r)) 



we obtain 



V{cP)<2{-l3o-S + k(3o)- 



\0) 



0, 



'(l + <^n(0))2 

since k^o = 5 + /Sq. 

The hypotheses of Theorem 1 are satisfied with a{u) = G{u), and it 
follows that xi = is stable in the case of this subsection. 



2.2 Stability properties of X2 

In this case, 

B = Po[n-{n-l)A]/A^ (17) 

where A = l3o{k - l)/6. 

As pointed out in Subsection 1.2, in this case the condition ^ must be 
fulfilled. 

We refine the discussion concerning the cases of stability given at the 
beginning of Section 2, for this concrete B. 

I. A. The condition B < and the definition of B imply n — {n — l)A < 0. 
This implies n > 1 and 

f(.-:)>^. ,18) 
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The condition B + 5 < leads to n — (n — k)A < 0, that impHes n > k 

and 

Obviously, the second inequality implies the first one. 

The sufficient condition of local stability is condition ()15p , with B given 
by (|17p . We remark that the condition |(5 + i?| < \kB\ is satisfied since it is 
equivalent to 6 + B > kB and this one is equivalent to ^(A; - 1) > 1, (the 
condition of positivity of X2). 

We have to study the behavior of the solutions at the extremities of the 
interval of stability. 

a) We consider the case 

^ ^ avccosiiS + B)/{kB)) ^ ^^^^ 

LOq 

This relation, together with 

ojocot{ujor) = -{6 + B), (21) 
and Wo > (from the definition of wq), imply 

L^O = V(fcS)2-(5 + S)2 (22) 

and that the pair fi* = 0, u* = represents a solution of ([9]). 

For later use we remark that, for B < 0, the relations (j2ip and ()22p 
(where, by the definition of ujq, ujQr € (0,7r)) together, imply relation ([20]) 
and again that the pair /x* = 0, w* = wq is a solution of ([9|). 

We assume that we vary one of the parameters, that we denote here 
by a, such that for a value a* the equality ([20]) is satisfied, and keep all 
other parameters fixed. We then obtain two complex conjugated branches of 
eigenvalues Ai^2(a) = fJ,{a)±iuj{a), such that Ai^2(a*) = ±«w*. If ^(a*) 7^ 
and the first Lyapunov coefficient of the reduced on the center manifold at 
a* equation is different from zero, then a Hopf bifurcation takes place in the 
center manifold. The sign of the first Lyapunov coefficient gives the stability 
properties of the solution at a* and of the periodic solution that occur by 
Hopf bifurcation. If the first Lyapunov coefficient is equal to zero, then a 
degenerated Hopf bifurcation takes place at a* . 

The construction of an approximation of the center manifold and the 
computation of the first Lyapunov coefficient (and thus of the normal form 
of the reduced equation) at a Hopf bifurcation point constitute the object 
of another paper of ours, [4j. 

b) The case r|(^ + i3| = 1, corresponds to the case ai = 1 from the paper 
of Hayes, [2]. In this case there always are eigenvalues with either positive 
or zero real part. The case of eigenvalues with zero real part (and all other 
with negative real part) corresponds to the case 02 = — 1 of [2j. By using 
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the relations between ai, 02 and p, q, (these are oi = —pr, 02 = qr) we find 
pr = —1, qr = —1, hence 5 + B = kB = —1/r. The first equality implies, as 
above, ^{k — 1) = 1 and it can not be satisfied in the zone of the parameters 
that we consider here. It follows that when r\S + B\ = 1, the solution X2 is 
unstable. 

Remark. Assume that we vary r and keep all the other parameters 
fixed. The conditions (|15p or (|16p for r are not as simple as they seem, 
because B is itself a function of r (being a function of k). Let us consider 
the function 

g{r) = T-\-{S + B{r))r) - arccos (^(^||^) • (23) 

If for a certain r* we have g{r*) = (that is the condition for the change 
of stability), in order to find whether a value ri in a neighborhood of r* is in 
the stability zone or not, we have to know the sign of g{ri), hence we have 
to study the monotony properties of function g in a neighborhood of r*. 

I. B. Since here B < 0, B + 6 > we must have 

^{k-l){n-k) <n. (24) 


The sufficient condition of local stability is condition (|16p , with B given by 

m- 

A point in the parameter space, satisfying 

arccos(( S + B)/kB) 
r = 

Wo 

is a Hopf bifurcation point, if, with the notations from I. A. a, ^(a*) 7^ 
and the first Lyapunov coefficient of the reduced on the center manifold at 
a* equation is different from zero. The stability of the solution in this case 
is given by the sign of the first Lyapunov coefficient. 

Remark, B < 0, 5 + B = 0, the solution is stable if and only if 
—kBr < Tr/2 while, for this case, the point kBr = —tt/2 is a Hopf bifurcation 
point. 

II. i? > 0. If n - 1 < 0, then i? > 0. If n - 1 > 0, then i? > is equivalent 

to 

f(.-l)<^. (2.) 

In this situation, all the eigenvalues have negative real part if and only 
if kB < 6 + B. This inequality is equivalent to 

^— B < 1 ^ - n - (n - 1^ <1 ^ A> 1 ^ ^ > 1, 

and this last inequality is already imposed (by the condition X2 > 0. Hence 
in the case i? > 0, X2 is stable. 
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3 Comments on the stability results in [51, [S] 



In order to compare our results with those of [5], [6], we define, for n > 1, 

r„ := ln<^ - - + 1 

and remark that r„ > ^^(5 < f^Q. Also for n > 1, relation (j'25|) implies 
that B > ^ r > Vn- This last condition is trivially accomplished if r„ < 
which is equivalent to > /3o- 

With these remarks we get the following situations for the sign of B. 

I. If n < 1 then B > 0. 

II. If n > 1 and > /3o then B > 0. 

III. If ra > 1 and < /3o then 

5 > for r„ < r < r^ax, 

B < for < r < r„ . 

This discussion allows us to follow the results of [5], [6] (the delay is 
there denoted by r). Those results have the following weak points. 

1. The results in [5] are presented for both equilibrium points simul- 
taneously, and this leads to imprecisions. As example, the affirmation at 
point (1) in [5j, pg. 238, is not true for xi = 0. Actually, the characteristic 
equation for this equilibrium point does not depend on n and thus for this 
point the condition n £ [0,1] is irrelevant. The condition of stability for this 
point does not depend on n. The ambiguity induced by using the plural 
"solutions" persists also at point (2) of [5J, pg. 238, leading to misunder- 
standings since the conclusions there can not refer to xi, as is seen from our 
Subsection 2.1. 

2. In the case B < 0, the sign of 5 + B is not considered in 0, pg. 238. 
To express the results of the analysis of the sign of 6 + B in terms of r is 
a little more difficult since inequalities ()19p and (|24p contain second degree 
terms in k. However the cases B + 5 > and B + S < are different in 
conclusions and they can not be eluded. 

3. The conclusions in [5], pg. 238, (2), b) seem to refer to the case 
B < 0, 5 + B > 0, but even for this case the result therein is not correct, 
since there the stability condition is 



r < 

instead of 

r < 



arccos{{B + 6)/ kB) 
^{kBY - {6 + W 

avccos{{B + 6) /kB) 



UJQ 

with coq defined in ([13]), as it is correct (condition (fTUj) ) 
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Since, in general, wq 7^ sj (kB)"^ — {5 + BY (equality holds, for i? < 0, 
only when relation (|20p holds), it is obvious that the domain of stability 
found in [5] is not correct (not even for the case B < 5 + B > Q). 

In [6j, pgs. 316-317, the results are basically the same as in [5], excepting 
the fact that it seems that the discussion refers only to X2 (but the plural 
"solutions" is used again). However, the observations from 2. and 3. above 
remain valid for [6j also. 
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